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Abstract 

In this paper we will introduce the methodology of analysis of the convex hull of 
the attractors of iterated functional systems (IFS) - compact fixed sets of self-similarity 
mapping: 

K= (J Ai{K)+U (1) 

i=l..n 

where Ai are some contracting, linear mappings. The method is based on a function 
which for a direction, gives width in that direction. We can write the self similarity 
equation in terms of this function, solve and analyze them. Using this function we can 
quickly check if the distance from K of a given x is smaller than a given distance or 
even compute analytically convex hull area and the length of its boundary. 



1 Introduction 



Fix X = W n with \\x\\ := Vx? x norm, A, : X — > X(i = l..n) are contracting matrices 

q := \\Ai\\ := sup < 1 

x: ||x|| = l 

and ti £ X for i = l..n are translations. 

The finite set of contracting mappings {x — > AiX + tj}i=i.. n is called iterated functional 
system. In [2] it is proved that it has an attractor: a unique compact nonempty set K, which 
is fixed point of I: 

I(K):= |J (MK) + U) (2) 

i=l. .n 

where A(K) := {Ax : x G K}. 



A precise analysis of this attractor is usually very difficult. It's common problem in 
e.g. computer graphics, image compression to find a good approximations from above of 
this set. There are known approaches to bound this set using spheres: Rice[l], Hart and 
DeFanti[3], Sharp and Edalat[Sj or boxes [Hj. 

In this paper will be introduced the methodology to bound it with convex set too, but 
this time - optimally: we will show how to construct its convex hull and how to use to get 
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quickly as precise approximation of the attractor as needed. 

We will reduce this problem to the solution of n — 1 dimensional functional equation for 
the width function, which can be easily approximated numerically and in some cases even 
found analytically. 

In Section 2 we will define the width function - it's a function which for a given direction, 
gives a position of orthogonal hyperplane bounding the set in this direction. This function 
completely defines the convex hull of closed set. 

To describe a convex set, it's better to use the radius function for given point, which gives 
for a given direction, length of segment from that point in that direction. We will show how 
to get the radius function from the width function. 

In Section 3 we will show how to change the self similarity equation for K to functional 

equation of its width function and that we can solve it numerically by a simple iteration. 

In Section 4 we will show how to use the width function in practice - for example to decide 

if the distance from the point x to K is smaller than a given number. 

In Section 5 we will show explicit formula for the width function for IFS with all Ai equal. 

Then we will concentrate on simple 2-dimensional IFS. We will show that it has a point of 

symmetry and using the width function - that its convex hull is built of triangles. So we can 

compute the area of the convex hull and the length of its boundary. 

Using the isodiametric inequality we can get an interesting trigonometric inequality: 

Vi< rG R, <pe[Q.M 52 I sin(j0)|r~' ? < -^7- 

j>0 



2 Width function 

Definition 1. For a bounded, nonempty set K C X we define width function h x : D — > R 
around x G X : 

h K ,x(d) := inf{/i : K C x + H(d, h)} 

where D := S m _i = {16 IR m : x T x =1} - directions, 
H(d, h) := {iel m : x T d < h} - halfplane. 

K will be usually fixed, so we will write h x = hx,x 
Of course h x (d) = h y (d) + (y — x) T d 
If xq is in the convex hull of K, h Xo is nonnegative. 

Obviously h is bounded and it is easy to show [1] that is continuous, and that it completely 
describes the convex hull of compact set: 

L = conv(A') = p| H(d, h K:Q (d)). (3) 
deD 

Fix a nonempty, compact, convex, set L C X and Xq G L. 
Definition 2. For L, xq as above we define radius function r : D — > M + U {0} around Xq: 

( r L,x (d) =) r(d) := sup{r : x + rd G L} (4) 
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Now we can analyze the relation between the radius and the width functions of L around 
the same point, say x = (h = h ). 
For any s > 



r(d) < s <^> V egD:d T e>0 sd G H(e, h(e)) O V eeD;d T e>0 sd e< /i(e) V egD:d T e>0 s < 



So 



r(d) 



inf 



Me) 



i&D:e T d>0 d T C 

If e fulfils infinium for d we say that e supports d. 



(5) 



We would like to use differential methods - we have to expand h for a moment in a 
neighborhood of the sphere. We would do it in the simplest way - take h(e) to ensure that 

e T S/h(e) = 

where x := x/\\x\\, V = (9 e< ) i=1 .. m . 

Assume now that /i(e) is differentiable for some e £ D. 
So the necessity condition for e G -D to support some d £ D from (EJ) is 

eKe<s> 3 AGR (d T e)V/i(e) - h(e)d = Xe 

d 1 e 



multiplying by e T , we get A = —h(e)(e T d) 

d = (d T e) 



s fVh(e) 



+ e 



Now we take d T e > to fulfil 



1, (ISI) becomes 



Vh{e) 



h(e) 



(6) 



3 Self similarity relation 

For any nonsingular matrix A : X — > X, translation t£X,x£X,d£D,a£ 

Ax + t G H(d, a) <^> x T A T d < a -t T d <^> x G H (A T d, a * f 

\ \\A "II 

Now look at the self-similarity equation (1): 

h (d) < a K c H(d, a) Vj A { K + UC H(d, a) V { K c H (^fd, 



(7) 



ta a - t J d \ a 

Pf4 ' 



^ V 4 /io(i4fd) < f=£j ^ max, ( d|| + tfd) < 



It's true for any a G M + , so we've get functional equation for the width function: 
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Observation 3. 



h (d) = max || Aj d\\h {A T d) + tjd. 



(8) 



In some cases we can solve this equation analytically, but numerical approximation should 
usually be enough. 

Consider C(D, M.) - the space of continuous functions D — > R with supremum norm: 
11/11 = swp deD f(d). 
Define I : C(D,R) -> C(D,M) 



I h (f)(d) ■.= max(\\Afd\\f(Ajd)+tfd 



Now for any f,g£ C(D, . 



\\h(f) ~ h(g)\\ < sup sup 

d&D r£[-R,R] 

< sup sup 

d&D r£[-R,R] 



max 



max 

i 



Afd\\(g(Afd)+r)+tJd 
max (||A^d||r) = i?max = ||/ — g\\c 



(9) 



\\Ajd\\g(Ajd) + tjd 



where R := \\f — g\\. 
c := max,- IIA- 



max,- c; < 1 . 



So because are contracting, Ih is contracting with c coefficient. 



< 



Using Banach contraction theorem, we get the unique fixed point of iteration (Q - the 
width function of our attractor. 

So to approximate the width function numerically, we can start from a constant function 
(the width function of a ball) and iterate (Q. 



4 Approximation of attractor 

In this section it will be shown how to use found width function to approximate K as precise 
as needed. 

To check if x G L (convex hull of K) we should check if 

x E L <-> \/ e <=D(x - xofe < h Xo (e) <-> \\x - x \\ < r Xo (x - x ) 
for some fixed Xq G L. 

We can check it immediately having the radius function, but finding this function from 
the width function ([6]), requires some smoothness - can be generally difficult, especially for 
numerical approximated functions. 

We will see that we won't loose much of precision, if instead of checking all directions for 
the width function, we will check only one: x — Xq 

L X0 := {x : (x - x ) T (x~^~Xo) < h Xo (x^Xo)} = {x : x T (x^~x ) < h (x^~x )} 

From fig. 1 we see that L XQ for a segment, where xq is one of its ending points, is a ball 
with the center in the middle of this segment. 
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Figure 1: L Xo for a segment (left) and some polygon. 

We can threat convex set(L) as the sum of such segments for all directions, so L xo is the 
sum of all such balls (in fact we can restrict to the supporting points of K). 
So we have rough approximations: 

L C L XQ c B(x , R xo ) 
p(L X0 ,L)<R X0 /2 

where p(A, B) := swp xeA mf ye B \\x — y\\, 

R xo : = sup deD h xo (d) (= inf r>0 {L C B(x ,r)}). 
If h Xo isn't constant (L is a ball) this approximation is better than with a ball. 

We could go closer to L this way, by checking more points, but we can use this additional 
checkings to came closer to the attractor - check if x G I k (L X0 ): 

x G I l+1 (L Xo ) ^3 iX e A(I\L X0 )) +ti^3i A~ l (x - U) G I\L X0 ) 

If we iterate this equivalence k times, we get an algorithm checking if x G I k (L X0 ): 

function near (x , k) 
{ int i=0; 

if (x T (x — Xo) > h (x — Xo)) return false else 
if(k=0) return true else 
{ do i++ until ((i>n) or (x E lm(Ai) and near(A^" 1 (a; — U) ,k-l))) ; 
if(i>n) return false else return true 

} 

> 

The condition with the image of Aj, is required for singular matrices and numerically 
cannot be fulfilled - we can just omit this matrices. 
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To evaluate the results of this algorithm, define 

C k = p(I k {L X0 ),K) 

Co < Rx a /V% is some constant, which should be approximated from above eg numerically. 
Let's make the iteration now: 

inf ||x - y\\ = inf inf \\x - U - A^y)]] < inf ||A|| mf \\Ar\x - *i) " v\\ (10) 

y£K i y&K i y&K 

So Ck < Ck-i maxj \\Ai\\ = Ck~i maxj q < Coc k . 

Now using (1101) we can alter this algorithm: if nearl(x,0 will return true, we are sure 
that x is nearer to K than a given distance I: 

function nearl(x, /) 
{ int i=0; 

if (x T (x — Xq) > ho(x — Xo)) return false else 
if (/ > Co) return true else 
{ do i++ until ((i>n) or nearl (Aj l (x — tA, l/ci) ) ; 
if(i>n) return false else return true 

} 

> 

Analogically we can construct for example an algorithm to answer questions like if 
I k (L)f)(I k (L)+t) = 0. 

5 Analytically solvable examples 

In this section we will show that in some cases we can compute width function analytically. 
Observation 4. In the case A\ = A 2 = ... = A n (= A) the width function realize 

h (d) = \\Ad\\h {Ad) + h*(d) (11) 

where h*(d) = maxitfd = h{ tl ,t 2 ,...,t n },o(d) 
Solution to this equation is: 

oo 

h (d) = ^\\A i d\\h*(A r d) (12) 

8=0 



Proof: pip is obvious from flb]) 
By induction over k, using (TTTT) : 



fe-i 



ho(d) = \\A k d\\h Q (A k D) + ll^ll^(^) 



i=0 
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A is contracting, h is bounded - the first term tends to 0. □ 

Now we will show an example how to get complete set of information about the convex 
hull L of K using the width function. 

Fix X = M 2 . We can threat X = C as the complex plane, multiplication by a complex 
number corresponds to a rotation and a scaling. 

We can identify the directions space D with angles: a = (cos a, sin a). 

Fix re i(t> = z E C : r = \z\ > 1, 2<n6N 
In the rest of this section we will analyze attractor, which can represent fractional part in a 
complex base system [8]: 

Kz= ( K + i ) ( 13 ) 

j=0..n-l 

Before investigating the width function, we observe that K has the symmetry point: 
Observation 5. Xq = |f5r ^ s the center of symmetry of K. 
Proof: For any point Xq, from f)13p 2x — K satisfies: 

z(2x - K) = 2x z - (K + i)= ^ 2x o z -n+l-K + k 

i=0.. n— 1 k=n— 1— i=0. .n— 1 

Hence if 2xoz — n + 1 = 2x then 2x$ = and we have: 

z(2x - K) = J2 ( 2x o " K ) + 1 

i=0..n-l 

Because of the uniqueness: 2xq — K = K, so x _ K = K — xq. □ 

We have h xo (a) = h(a) = h(—a) = (h (a) + h (—a))/2 
Now h*(a) + h*(-a) = (n - l)r _1 | cos(a + 0)|, gives 

2h(a) = {n - 1) ^ r~ j \ cos(a + j<f>) | (14) 
i>o 

Now assume that = 7r^ (I, k E N) 
2h(a) = (n - 1) I cos (« + 34>) I E r " i "* < = ^ r "' 1 C ° s(a + J0) 1 (15) 

j'=l..fe i>0 j=l..k 

In this case, our width function is differentiable everywhere except finite number of 
points. We could use (jSJ) (Vh(a) = (— sin a, cos a)h'(a)) for differentiable points of h 
around indifferentiable, but instead we will show that this points corresponds to straight 
lines on the boundary of L (convex hull of K). In fact, it will occur that L is a polygon. 
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Figure 2: Analysis of the behavior of width function around a segment on the boundary of 
a convex set (left) and construction of the convex hull in case: z — i + 1, n = 2 (right). 

Look at the triangle on fig. 2 - width function will have indifferentiable minimum in - 
in a neighborhood of 0: 

, / \ _ / dcos(j) > a < e 
~ I dcos(5) > a > -e 

h'(0 + ) = rfsin(a + 7 ) = b , /i'(0") = -esin(/3 + 8) = -c (16) 

So we see: 

1. If a vertex in some direction 7 is supported by some range of directions - in that range 
h(a) = g?cos(7 — a). Thanks to the uniqueness of h, implication above is equivalence. 

2. In a direction in which the edge of the convex hull contains straight segment, iff the 
width function have local minimum and is indifferentiable. 

We can construct this segment using (fIBI) . 

Now look at (TloT) - there is a finite number of indifferentiability points and between them 
function is in the form of h(a) = J2i=o k-i a * cos ( a + h)- Expanding cosinus of sum, we see 
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that it can be written in the form h(a) = a cos(a + b) for some a, b e R. So it supports on 
one point - common point of succeeding triangles - our convex hull is polygon. 



We will construct it now - call T(a, b, c) - triangle as in fig. 2 
This triangles should be rotated by angles of indifferentiability (fij := ir/2 — j<p for j = \..k 

dj = h(4>j) 

h + <h = h'(<pp - h'(<j>j) = l -{n - 1)^-^( COS '(0+) - cos'(07)) = („ - l)y^ 
b .- Cj = h 'ty+) + h 'tyj) = J!L^ r-^|cos(a + i0)| a=0 , 

i=X..k,i^=j 

L = x + |J e^iTia^b^c^U-Tia^bj,^)) (17) 
j=i..k 

Notice that if we wouldn't gather elements with the same value of cosinus in (|T5|) we can 
look on the construction above as(now j = I..00): 

aj = h{<f>j) 

bj + cj = ti{<ft) - ti((j>j) = hn - l)r- J '(cos'(0+) - cos'(0-)) = (n - l)r" J 
bj ~ Cj = h'((pj) + h'((j)j) = (n - 1) ^2 r_l ^"l c °s(a + «0)U=^ 

i>0,i^j 

L = x + |J (T(a„ 6,-, Cj ) U -T(a„ 6,-, c,)) (18) 

j>0 

In form ([TBI) triangles from (TTTj) are constructed from infinity many triangles with disjoint 
interiors. 

We will use this form to construct L in <fi ^ 7rQ case. 

Take (ff G 7rQ, lim^oo = series. We construct polygon for each <p l . It's easy to check: 

<t>)^<t>j u) ■ % b) ■ b, r) > Cj 

for any j > 0. 

So 018p gives convex hull in this case. 

Its boundary contains countable number of segments. 

We can now find a formula for length of the boundary(5) of L and it's area(A): 

n — \ 



r — I 



B = 2(n- l)^2r~ j 

j>0 

A = -(n- l) 2 ^V 4 ^r^'|cos((z-j> + 7r/2)| = -(n - 1) 2 | sin(0)|(r- 2 + r" 4 + ...) + 

j>0 j>0 

+ {n-l) 2 I sm( V( f))\(r~ 2 - v + r~ 4 ~ v + ...) = ^ | sm(v(/>)\r- v 



v.=\i—j\=l..oo v>0 
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It's interesting that B doesn't depends on 0. 



We can get an interesting trigonometric inequality from this formulas, namely we know 
(eg [7]) that for given length of boundary, the largest area has circle, so A < B 2 /Ait: 

Observation 6. Ylj>o I sin(j0)|r"- ? < ^^zj f or an V ?" > 1, G [0, 2ir]. 

Digression: to analyze points of indifferentiability in higher dimension - they will cor- 
respond to pyramids, which base can be analyzed by cutting space with two-dimensional 
planes - jumps of derivatives in different direction gives the width function for the base of 
pyramid. As above, we can also do it in different way - by analyzing differentiate points 
around indifferentiable one. 



6 Conclusion 

We have shown that self similarity equation can be written in terms of function describing 
some of its properties, like width in any direction. Other property which can be written in 
that way and can give some interesting results can be f(x) := fi(K fl (K + x)). 
Obtained functional equation can be usually approximated numerically and used to approx- 
imate, analyze our set. 

The width function can be used for example to exclude some points or to find convex hull 
of it and calculate some of its properties. 
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